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SUMMARY 


Finite-difference methods for unsteady transonic flow frequently use simplified 
equations in which certain of the time-dependent terms are omitted from the governing 
equations. This paper derives kernel functions for two-dimensional subsonic flow 
which provide accurate solutions of the linearized potential equation with the same 
time-dependent terms omitted. These solutions make possible a direct evaluation of 
the finite-difference codes for the linear problem. Calculations with two of these 
low-frequency kernel functions verify the accuracy of the LTRAN2 and HYTRAN2 finite- 
difference codes. Comparisons of the low-frequency kernel-function results with the 
Possio-kernel-function solution of the complete linear equations indicate the 
adequacy of the HYTRAN approximation for frequencies in the range of interest for 
flutter calculations. 


INTRODUCTION 

Time-dependent finite-difference methods are receiving increased use for aero- 
elastic calculations in transonic flow. These methods are implemented by computer 
codes that solve various approximations to the time-dependent nonlinear potential 
equation. For example, the widely used LTRAN2 code (ref. 1) provides a solution to 
the two-dimensional small disturbance potential equation in which the only time- 
derivative term retained is the first derivative term in the differential equation. 

A more accurate approximation is made with the HYTRAN2 code (called LTRAN2-NRL in 
ref. 2), which also includes the time-derivative terms in the boundary conditions and 
pressure calculation. In the past, the adequacy of these approximations has been 
assessed by comparing results from these codes with calculations based on an analytic 
solution of the complete linear equations. In contrast with the finite-difference 
codes, these linear calculations have included the second time derivative in the 
differential equation. 

In order to more easily assess the significance of omitting some of the time- 
dependent terms from the equations and to evaluate the adequacy of the finite- 
difference codes in solving these equations, kernel-function solutions to the simpli- 
fied equations have been derived. This paper gives the derivation and compares cal- 
culations using these new kernel functions with classical results from the Possio 
integral equation. In addition, the accuracy of the finite-difference codes in solv- 
ing the linearized simplified equations is demonstrated by comparing them with the 
kernel-function solutions. This comparison verifies the linear theory aspects of the 
code and may be used to assess the adequacy of the computational grid. 


SYMBOLS 

b airfoil semichord, m 

Cp pressure coefficient 

ACp lifting pressure coefficient 


slope of lift coefficient 



slope of pitching-moment coefficient about pitch axis 


Cm„ 

Hankel function of second kind for order i 
K kernel function 

k reduced frequency, bto/U 

M free-stream Mach number 

N integral (see eg. (A23) or (A43)) 

P pressure term in integral equation, ACp 

P^V2,-1/2) Jacobi polynomial 

r scaled x-variable, M2k(x - 

s scaled y-variable, M^ky/p^ 

t dimensionless time, Ut'/b 

t' time, s 

U free-stream velocity, m/s 

u scaled x-variable, k(x - C)/P^ 

V scaled y-variable, ky/p^ 

w dimensionless downwash, w'/U 

w' downwash velocity, m/s 

X dimensionless x-variable, x'/b 

x' streamwise coordinate, positive downstream from midchord, m 

y dimensionless y-variable, Py'/b 

y' coordinate normal to free stream, m 

z dimensionless mode shape, z'/b 

z' airfoil mode shape, m 

p = ''jl - 

9 angle variable (see sketch B) 

X, n, T integration variables 

V parameter (see eq. (A10)) 
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dummy x-variable 


? 

p radial variable (see sketch B) 

a Fourier transform variable 

(|) dimensionless velocity potential, <|)*/bU 

({)' velocity potential, m /s 

cb dimensionless acceleration potential, - 4 C 

2 p 

0 ) oscillation frequency, rad/s 

Subscripts: 

H HYTRAN 

L LTRAN 

P Possio 


ANALYSIS 

In this section the governing differential equation, the boundary conditions, 
and the integral equation are presented for the classical Possio formulation and for 
the HYTRAN and LTRAN approximations to it. The labels HYTRAN and LTRAN are used for 
the kernel-function formulations which contain the same simplifications of the equa- 
tions as are used in the HYTRAN2 and LTRAN2 finite-difference codes. The details of 
the derivation are given in the appendix. 


Differential Equations and Boundary Conditions 

The problem is formulated in terms of the linearized, small-disturbance, pertur- 
bation velocity potential. The differential equations and boundary conditions are 
derived in section 5-1 of reference 3. They are given here in terms of dimensional 
variables. 


Complete linear equations .- In dimensional variables, indicated by primes 
throughout, the perturbation potential satisfies the linear differential equation 


(1 - 




Vy 



M 


U 




= 0 


( 1 ) 


where M and U are free-stream Mach number and velocity, respectively. The pres- 
sure coefficient is 


C 

P 



+ <!>;.) 


(2) 
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For the small perturbations treated herein, the lifting part of the problem is anti- 
symmetric in y' and we can limit our attention to the upper-half plane. For an 
airfoil motion given by z'(x',t') = 0, the airfoil flow-tangency condition is 

w' (x' ,0''’,t' ) = Uz^, + zj_, (3) 

where w' is the downwash velocity on the airfoil and b is the airfoil semichord. 
The remaining boundary condition on the x*-axis is to require zero pressure discon- 
tinuity; that is 

AC (x’,0,t') = C (x’,0“,t’) - C (x',0'^,t’) = 0 fix'! > b) (4) 

P P P V I I y 

In addition, the Kutta condition must be imposed at the trailing edge to insure a 
unique solution; this is enforced by requiring 

tC (b,0,t') = 0 (5) 

P 

Finally, for this time-dependent problem, the Sommerfeld radiation condition requires 
outgoing disturbances at infinity. (See sec. 6-4 of ref. 3.) 

HYTRAN equations .- The HYTRAN formulation differs from that for the linear equa- 
tions in the preceding paragraph only in that the 4>t't' term is omitted from the 
differential equation; that is, for this case, one solves 


(1 - 


M^) 


V'y' 


2M 

U 


^x't* 


= 0 


( 6 ) 


The boundary conditions are unchanged. 

LTRAN equations .- The LTRAN formulation solves the same differential equation as 
does the HYTRAN formulation (eq. (6)). The boundary conditions are further simpli- 
fied by omitting the time -dependent terms; that is, the pressure coefficient is given 
by 



and the airfoil boundary condition is 


(7) 


w' (x' ,0'^,t' ) = Uz^, 


(|x* I < b) (8) 
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Integral Equations 


The equations given in the preceding sections are formulated in the appendix by 
use of integral transforms of the acceleration potential forms of the equations. In 
each case, the solution is obtained from an integral equation relating the known 
downwash on the airfoil to the unknown lifting pressure, namely 


w(x) 



P(5) K(x-5) d^ 


( 9 ) 


where nondimensional variables are now used, and w is the downwash, P is the 
lifting pressure, and K is the kernel function. The time dependence is removed by 
means of the Fourier transform and the resulting equations are for harmonic motions. 

The three formulations differ only in the form of the kernel function. The 
forms given in this paragraph are in terms of singular integrals which must be inter- 
preted as finite-part integrals in Hadamard's sense. Computational forms are given 
in the appendix. The results, expressed as integrals of the Hankel function, are as 
follows: 

Possio kernel function; 


K 


P 


iMk -ik(x-^) 


J 


k(x-5)/p" 


e 


h"’(m|x|) 


dX 

w 


( 10 ) 


HYTRAN kernel function: 


K 


H 


iM^k -ik(x-F) 

8p ® 


/ 


k(x-^)/p" 


iX 

e 



(11) 


LTRAN kernel function: 


K. = - 


. 2 
iM k 

8p 


I 


M^k(x-C)/p^ 


.,( 2 ) 


H 


1 


(IM) 


dX 

W 


( 12 ) 


Note that the HYTRAN kernel function may be obtained from the Possio kernel func- 
tion by replacing M with in the two places M appears explicitly (not in 


p = \ll - M^). 
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Finally, equation (12) demonstrates a transonic similarity law for the LTRAN 
formulation. One observes that K^/P depends on a single parameter M^k/p^ and the 
variable x - The integral equation (9), with a slight rewriting is 


w(x) 



ppf 


(13) 


so that for this equation pp also depends on the single parameter M^k/p^. It 
follows that the aerodynamic forces obey the same similarity law. 


Solution of Integral Equation 

Integral equation (9) was solved by a collocation procedure. The unknown 
pressure was approximated by a sum of 64 weighted polynomials 


P(!i) 


64 

I 1 - ^ sin [(n - 1/2) arccos 

II + Z 2^. ^n sin [(1/2) arccos 

n=1 


(14) 


with unknown pressure coefficients aj^. The ratio of the sine functions is pro- 

( 1/2 - 1 / 2 ) 

portional to the Jacobi polynomial ^ ' (ref. 4, ch. 22) appropriate to 

the pressure weight function ^ ( 1 - ^)/( 1 + ^) . (This weight function enforces the 
Kutta condition.) Equation (14) is substituted into equation (9), which leads to 


w(x) 



1 - g; sin [(n - 1/2) arccos g;] 
1 + ^ sin [(1/2) arccos 


K(x-^) d^ 


(15) 


The integral is evaluated by the Jacobi-Gauss quadrature appropriate to the weight 
function for a set of 64 collocation points at which the downwash is specified. This 
procedure leads to a set of 64 simultaneous equations which are solved for the 
unknown pressure coefficients a^^. The choice of 64 collocation points gives solu- 
tions which have converged to plotting accuracy over the range of frequencies and 
Mach numbers considered. 


RESULTS AND DISCUSSION 

Calculations of aerodynamic forces for an airfoil oscillating in pitch were made 
with each of the three kernel functions. Comparisons among these kernel-function 
results are presented which show the influence of Mach number and frequency on the 
lift and pitching moment. A second comparison is made of calculations from the 
kernel functions with those from two finite-difference codes. 
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Kernel-Function Results 


Calculations were made for a flat-plate airfoil oscillating in pitch about the 
quarter-chord axis at Mach numbers 0.7, 0.8, 0.9, and 0.95. The slope of the lift 

coefficient c (fig. 1) and pitching-moment coefficient c (fig. 2) are shown 

a 

as functions of reduced frequency for each of the three kernel functions. In each 
figure, both the real part (in phase with the motion) and imaginary part (in quadra- 
ture with the motion) are shown. 

In general, the agreement between calculations with the Possio and HYTRAN kernel 
functions is good, with the real parts being in better agreement than the imaginary 
parts and the lift being in better agreement than the pitching moment. In particu- 
lar, for frequencies in the range of interest for flutter (k < 0.2, say), the two 
results are quite close together. In contrast, the LTRAN-kernel-function results 
differ greatly from those of the other two kernel functions even at quite low reduced 
frequencies. (An exception occurs in fig. 2(a) for which the curves for the real 
part cross at k = 0.58.) For example, at M = 0.8, the imaginary part of the LTRAN 
c-^^ (fig. 1(b)) differs by over 70 percent from the Possio value when k = 0.1. 

Flutter calculations using each of the three kernel functions are reported in 
reference 5. Those calculations show the agreement between Possio- and HYTRAN- 
kernel-function flutter results and the disagreement between these and the LTRAN- 
kernel-function results that might be expected in light of the comments of the 
preceding paragraph. 


Comparison With Finite-Difference Calculations 

The finite-difference codes of references 1 and 2 have been widely used for the 
unsteady flow analysis of airfoils at transonic speeds. These codes provide numeri- 
cal solutions to approximations of the nonlinear, small disturbance potential equa- 
tion for time-dependent transonic flow. The codes, LTRAN2 and HYTRAN2, omit the same 
time- dependent terms from the equations as do the corresponding LTRAN and HYTRAN 
kernel functions, as described in the section "Analysis." 

These inviscid, finite-difference codes have been evaluated in the past in two 
ways: by comparison with experiment or by comparison with exact linear theory. The 

difficulty with the former includes, among other things, the effects of viscosity and 
wind-tunnel walls. The problem with the latter comparison is that the equations used 
for the linear theory corresponded to the Possio kernel function. The development of 
the LTRAN and HYTRAN kernel functions herein has overcome this difficulty and has 
made a comparison of the finite-difference and analytical solution of the same equa- 
tions possible. 

An example of such a comparison between finite-difference calculations and the 

present analytic results is given in figure 3. The lift coefficient c is shown 

a 

as a function of reduced frequency for an airfoil oscillating in pitch about a mid- 
chord axis at M = 0.8. The curves give the analytic results for each of the three 
kernel functions. The symbols show the solutions obtained from the finite-difference 
code of reference 2. This code permits the use of either the LTRAN2 or HYTRAN2 for- 
mulation of the problem. For the present comparisons, these codes were run with the 
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nonlinear terras omitted for a flat-plate airfoil to correspond with the kernel- 
function results. These solutions for harmonic oscillation were obtained by running 
the time-dependent codes through the three cycles of oscillation required for the 
transients to die out. There were 360 time steps per cycle, and a Fourier transform 
of the last cycle of oscillation was used to obtain the harmonic loads. The original 
79 X 99 LTRAN2 grid was used with 33 points on the airfoil. 

The agreement between the finite-difference and analytic solutions shown in fig- 
ure 3 is very satisfactory. This comparison demonstrates, at least for this linear 
case, the accuracy of the algorithm, correctness of the coding, and adequacy of the 
boundary condition treatment, mesh spacing and extent, and time step size for the 
finite-difference methods. 


CONCLUSIONS 

The kernel functions for two levels of approximation to the linear equations of 
subsonic compressible flow have been derived. These functions were used in an 
integral-equation solution and compared with results from the classical Possio inte- 
gral equation in order to assess the limitations of the approximations. In addition, 
comparisons were made with two widely used finite-difference solutions of the approx- 
imate equations. 

Specific conclusions include 

( 1 ) The approximation which retains the time-derivative terms in the boundary 

conditions ( HYTRAN-like) agrees well with the complete solution over the 
range of frequencies of interest in flutter calculations. 

(2) The approximation which omits the time-derivative terms in the boundary 

conditions (LTRAN-like) is adequate only at very low frequencies. 

(3) The LTRAN2 and HYTRAN2 time- dependent finite-difference codes yield accurate 

solutions of the corresponding approximate equations, at least for the 
linear case. 

The analytical solutions derived herein provide objective standards for assess- 
ing attempts to enhance the efficiency of the finite-difference codes by techniques 
such as reducing mesh extent and increasing time step size. It is possible, of 
course, that the conclusions reached here may need to be modified for the nonlinear 
transonic-flow applications for which the finite-difference codes were developed. 
However, the kernel-function results can still be used to determine limits for the 
finite-difference parameters which may not safely be exceeded. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 2, 1982 
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APPENDIX 


SOLUTION OF DIFFERENTIAL EQUATION 

The derivation of the integral equation from the differential equation is pre- 
sented in this appendix. The velocity-potential statement of the problem is refor- 
mulated in terms of the acceleration potential and solved by integral transform 
methods. The kernel functions for each of three levels of approximation are obtained 
and forms useful for calculation are derived. 


Complete Linear Equations (Possio Kernel Function) 

Kernel-function derivation .- The derivation of the Possio kernel function uses 
the acceleration potential to solve the linear differential equation. A derivation 
by superposition of singularities is given in section 6-4 of reference 3. The nota- 
tion used herein is similar. The complete linear equations (eqs. (1) to (5)) are 
expressed in terms of the acceleration potential (j; and nondimensional coordinates. 
Let 


X' = bx '' 

y* = by/p 
t' = bt/U 
z' = bz 
w = w'/U 
(})' = bU(l) ^ 




(A1) 


and define the dimensionless acceleration potential to be 


{b = — !-c = <b + <b 

^ 2 p ^x ^t 


(A2) 


This expression may be inverted to express the downwash in terms of the potential as 


w(x,y,t) = P(|)^ ^ f d\ 


(A3) 
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which can be verified by differentiation. On the airfoil, when the antisymmetry of 
the pressure is taken into account. 


P(x,t) = ACp = Cp(x,0”,t) - Cp(x,0'^,t) 


= -2Cp(x,0'^,t) 


(A4) 


The time dependence is removed by performing a Fourier transform on t. The result 
is formally equivalent to specifying simple harmonic motion by letting 


(Kx,y,t) = (Mx,y) e^’^^ 


(A5) 


and similarly for other dependent variables, where 


k = ba/U 


(A6) 


In dealing with some of the integrals in the following analysis, it is necessary to 
apply analytic continuation arguments in cases for which certain integrals do not 
exist for real values of k. In particular, it is convenient to suppose that k has 
a negative imaginary part. 

Collecting these results, noting that ([j and satisfy the same differential 
equation, yields the following problem to be solved for y > 0; 


{|; + cl; - k(M/p)^(i2cl; - kcl,) = 0 

XX yy ^ X 




cKx,0'^) 



(|x| < 1)1 
(|x| > l)f 


(A7) 


w(x,y) 


p r (X,y) d\ 

oo ^ 




The problem is treated by using the Fourier transform 


(M o,y) 



-ic3x 

e 


(};(x,y) dx 


(A8) 
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The differential equation becomes 


dj - = 0 

yy 


(A9) 


in which 


- k(M/p)^(2a + k) 


(A10) 


The solution is (the branch being chosen for which Re(v) > 0) 


o , y ) = e CT,0^) 


1 - w 

= 4 ® 


r ' e-"^ 

-'-I 


P(^) dl 


(A11) 


The Fourier inversion formula 


‘h f 

C 


e^^*^ ^(a,y) da 


(A12) 


is used with equation (All) and the last of equations (A7) to obtain the integral 
equation 


w(x,y) 



P(^) K(x-C,y) dl 


(A13) 


in which the kernel function is 


K(x-C,y) 


= 1 - 1 - 
8x 9y 



-ik(x-X) 

e 


r” i(X-C)o -vy 
I e ^ e da d\ 


(A14) 


The integral-equation formulation of the original problem (eqs. 
setting y = 0 and using the boundary condition (see eq. (3)) 


(A7)) is obtained by 


w(x,0"'') 


Zj^(x) + z^(x) 


(A15) 
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in which z(x) is the nondimensional mode shape. The inner integral in equa- 
tion (A14) is evaluated as follows. Let 


2t 






a = a + ) 

P 

ff + 

X = 5 + -t- 


P^V 

^ k ^ 


(A16) 


Equation (A14) becomes (dropping the bars for convenience) 




-ip^u ^ iX C iXa -v^a^-M^ 
&v J ® I ® ® 

•'■■00 •' — CO 


do dX 


(A17) 


where now yy = v\[a^ - M^. Taking note of the parity of the inner integral yields 




■ip\ ^ r ^ix f " 

a' 1, ° Jo * 


cos Xa da dX 


(A18) 


The inner integral is a Fourier cosine transform which may be evaluated by using 
equation (1.4.26) of reference 6 and equations (9.6.27) and (9.6.4) of reference 4 
to obtain 


K(u,v) 


ik -ip^u 
8p ® 





(A19) 


The choice of the Hankel function of the second kind is consistent with the 
Sommerfeld radiation condition. For v > 0, the integral is well defined. To obtain 
the desired kernel function, one must let v ->• 0 . Using equations (9.6.1) and 
(9.6.27) of reference 4 yields 


lim 5^ 

v->0 . 2 
ov 




(A20) 
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This result is substituted into equation (A19) to give the final expression for the 
Possio kernel function 


Kp(u) 


iMk -iB^u 
8p ® 



i\ 

e 




dX 

W 


(A21) 


which is equation (10). 

This form for the kernel function is not useful for computation in that the 
integrand has a second-order pole at the origin. It is possible to derive the form 
given in equation (6-112) of reference 3 from it, but the form in the reference also 
involves an integral of the Hankel function. In the following section, a form is 
derived which involves only a standard quadrature of an elementary function. 

Computational form of kernel function .- The integral in equation (A21) does not 
exist for u > 0. This problem arose in taking the limit v = 0. Nevertheless, a 
useful computational form may be derived from it by using the integral definition of 
the Hankel function and deforming the path of integration in the complex plane. The 
formal analysis can be made rigorous by remembering that k may be assumed to have a 
negative imaginary part. 

Equation (A21) may be rewritten as 


Kp(u) 



-ip^u 


N(u) 


in which 


N(u) 



iX 

e 




dX 

w 


(A22) 


(A23) 


It is convenient to treat positive and negative u separately by using the sub- 
script + for u > 0 or - for u < 0. For u < 0, let X = -X and 


N = N 



-iX 

e 




For u > 0, 


(A24) 


N = N + N 
00 + 


where 


N 


+ 


f" iX „(2) . . , dX 
J e (MX) - 

*'|u| 


(A25) 


(A26) 
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and 


N 

00 




d\ 

w 


= 2 


/: 


(2) 

h1| ^ (MX) 


cos X - — 
X 


The integral in equation (A27) obviously does not exist. However, remembering 
source of the difficulty, one can use equation (A20) to rewrite it as 

N = - ^ ^ f hJ^^^(m\|x^+v^) cos X dX 

M v^ g^2 Jo 0 V / 

= -iZP 

M 


by equation (9.6.4) of reference 4 followed by equation (1.13.43) of reference 
Writing equations (A24) and (A26) compactly gives 


N - f r (MX) 

* J|u| ' ^ 

From equation (9.1.6) of reference 4, 


Using this result with equation (9.1.24) of reference 4, the integral becomes 


N 


± 



^i(±1-Mp) X 


dX d^ 


2M iu' f 

J, 


^-iM|u|fi N(|I^ - 1) 

1 + Mp 


d[i 


The integrand has poles at ^ = ±1/M and branch points at n = ±1. The path 
integration is deformed in the complex plane as indicated in sketch A. 


(A27) 

the 

(A28) 

6 . 

(A29) 


(A30) 

of 
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The integral along the arc vanishes since Re(iM|u||j.) > 0 there. The integral along 
the real axis can therefore be replaced by an integral parallel to the imaginary 
axis. Let p. = 1 - iO/ then 


a 


X 

m|u 


and 


N 


± 


i2 i(1+M)u 

e 

nMu 



-X ± i2Mu 
T - i( 1 + M)u 


(A31) 


Substituting this result and equations (A24) and (A25) into equation (A22) yields the 
computational form of the kernel function 


Kp(u) 


-k e 


2 

-iP u 



2 

1 - M 1 i(1+M)u 

SP 4Tipu ® 



fc e"' 


'jr ± i2Mu 
X - i( 1 + M)u 



(A32) 


The integral may be evaluated as a Laguerre-Gauss quadrature with weight function 
^ e"\ 
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HYTRAN Equations 

The HYTRAN formulation uses the same boundary conditions as the Possio formula- 
tion. However, the differential equation (eq. (1)) omits the <l>t't' term for the 
HYTRAN formulation. In the dimensionless acceleration potential variable of this 
appendix, the differential equation is 


(Jj + - i2k(M/p)^(P = 0 (A33) 

XX yy X 

for the case of harmonic motion. The derivation can be carried through in a parallel 
fashion to that given for the Possio kernel function. One finds 




2 2 
iM k “i6 u 

sp * 


f" h',2>(„2|M) 

•'—00 


d\ 

w 


(A34) 


which is equation (11). As noted in the text, the HYTRAN kernel function differs 
from the Possio kernel function in that the HYTRAN equation has in place of the 

M that appears in the Possio equation. 

The computational form of the HYTRAN kernel function is 


K„(u) = 


-k e 


2 

-ip u 


u \ )ji - -4 . . , .-..2, 


M , 1 i(1+M )u 

^ -Rj "8P ' ^ 


I 




-T IT ± i2M u 


2 

0 T-i(1+M)u 


dT 


(A35) 


LTRAN Equations 

The only time derivative in the LTRAN equations occurs in the differential equa- 
tion. The differential equation is thus the same as the HYTRAN equation, but the 
boundary conditions and pressure coefficient use the quasi-steady expressions. 

Kernel-function derivation .- The acceleration potential formulation is used 
again with the solution effected by use of the Fourier transform. The problem state- 
ment (eqs. (A7)) becomes 


(|; + c[j - i2k(M/B) 4> = 0 

XX yy ^ X 


+ fi 

(Kx,o^) = \ ^ P 


w(x,y) = p 


— 0 


(|x| < 1) I 
(l-l > 1) f 


(A36) 
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where now the acceleration potential is 


<Hx,y) = 


- C = ([) 

2 p X 


(A37) 


The most significant difference between the HYTRAN and LTRAN formulations lies in 
the absence of the convective exponential term from the integral in equations (A36) 
(cf . eqs. (A7) ) . 

The formal integral transform solution is carried out in a manner similar to 
that for the Possio formulation. The result is 


K (u) 

Li 


iM^k 

8p 



iX ,,(2) 


(|xl) 


d\ 

w 


(A38) 


which is equation (12). 

Computational form of kernel function .- In contrast to the Possio and HYTRAN 
kernel functions, the LTRAN kernel function can be expressed in terms of Hankel 
functions without an integration. For this case, the analog of equation (A17) is 






Let 


K = 

2 - 

a = M a 


u = r/M 


V = s/M 


and dropping the bars, one obtains 


K(r,s) = 


2 , 

M k 

Sup 


•'—CO •' — CO 


.( 0+1 )\ 


dX do 


(A39) 


(A40) 


(A41) 


Performing the inner integration, taking the indicated derivative, and setting 
s = 0 (y = 0) yields 


K(r) 


2 

iM k 
8n:P 



ire 

e 


0 - 1 
0 + 1 


d0 


(A42) 
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The branch of ^o^-1 has been chosen for which the real part of the radical is 
positive as |Re(a)| •> “. 

The integral (eq. (A42)) is evaluated by deforming the path in the complex 
plane. The branch chosen is defined in sketch B. 


a-plane 



For -TT < 


2 ^(01 + 60 V 2 

0 <_ 6^ < 2TT, \ a - 1 = 'VP 1 P 2 ^ ^ 

1- 5 l — n ' i(6. - 6 )/2 

-tr < e < TT, J-E — L = J — L e 1 

2- 1(0 + 1 Kp_ 


Sketch B 


It is convenient to treat positive and negative r separately. Let 


ira , a - 1 , 
e \ — r do 

'0+1 


For r > 0 , 


Re (ira) < 0 


in the upper half-plane. The path of integration may therefore be defoirmed into the 
path shown in sketch C. The integral around the branch point vanishes since the 
integrand is of order there. Below the cut. 
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and above the cut. 


0 


1 


2 


0 


The integrals above and below the cuts combine to give 


N^(r) 


= -2 



-ira 

e 


a + 1 
a - 1 


da 


(A44) 


For r < 0 , 


Re (ira) < 0 


in the lower half -plane and the path of integration may be deformed as shown in 
sketch D. Again the integral around the branch point vanishes since the integrand is 
of order \[^ there. Below the cut. 


and above the cut 
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One obtains 


N ( T) 



(A45) 


The two results for 



(eqs. (A44) and (A45)) may be written compactly as 


N (r> - 

± ^ a + 1 

- .2 f\--\r\a 1 ± g 

•*1 


do 


•= -21 ± 




(A46) 


The integral has been cast in the form of an integral representation of the Hankel 
function (eq. (9.1.24) of ref. 4). Thus, 


N 


± 


iitfl ± i 



Hj^’drl) 


(A47) 


Taking the derivative (eq. (9.1.28) of ref. 4) and substituting the result in equa- 
tion (A42) yields the final expression for the LTRAN kernel function: 


Kj^(r) - 


2 

M k ir 
8p ® 




where 


r 


M^U 


M k(x - 


P 


i) 


(A48) 


This form for the LTRAN kernel function is very convenient since Bessel functions are 
relatively easy to compute. 
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